3.1.99 \(\int \frac {a+b x^2}{(1+x^2+x^4)^2} \, dx\) [99]

3.1.99.1 Optimal result
3.1.99.2 Mathematica [C] (verified)
3.1.99.3 Rubi [A] (verified)
3.1.99.4 Maple [A] (verified)
3.1.99.5 Fricas [A] (verification not implemented)
3.1.99.6 Sympy [C] (verification not implemented)
3.1.99.7 Maxima [A] (verification not implemented)
3.1.99.8 Giac [A] (verification not implemented)
3.1.99.9 Mupad [B] (verification not implemented)

3.1.99.1 Optimal result

Integrand size = 18, antiderivative size = 119 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=\frac {x \left (a+b-(a-2 b) x^2\right )}{6 \left (1+x^2+x^4\right )}-\frac {(4 a+b) \arctan \left (\frac {1-2 x}{\sqrt {3}}\right )}{12 \sqrt {3}}+\frac {(4 a+b) \arctan \left (\frac {1+2 x}{\sqrt {3}}\right )}{12 \sqrt {3}}-\frac {1}{8} (2 a-b) \log \left (1-x+x^2\right )+\frac {1}{8} (2 a-b) \log \left (1+x+x^2\right ) \]

output
1/6*x*(a+b-(a-2*b)*x^2)/(x^4+x^2+1)-1/8*(2*a-b)*ln(x^2-x+1)+1/8*(2*a-b)*ln 
(x^2+x+1)-1/36*(4*a+b)*arctan(1/3*(1-2*x)*3^(1/2))*3^(1/2)+1/36*(4*a+b)*ar 
ctan(1/3*(1+2*x)*3^(1/2))*3^(1/2)
 
3.1.99.2 Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 0.16 (sec) , antiderivative size = 147, normalized size of antiderivative = 1.24 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=\frac {x \left (a+b-a x^2+2 b x^2\right )}{6 \left (1+x^2+x^4\right )}-\frac {\left (\left (-11 i+\sqrt {3}\right ) a-2 \left (-2 i+\sqrt {3}\right ) b\right ) \arctan \left (\frac {1}{2} \left (-i+\sqrt {3}\right ) x\right )}{6 \sqrt {6+6 i \sqrt {3}}}-\frac {\left (\left (11 i+\sqrt {3}\right ) a-2 \left (2 i+\sqrt {3}\right ) b\right ) \arctan \left (\frac {1}{2} \left (i+\sqrt {3}\right ) x\right )}{6 \sqrt {6-6 i \sqrt {3}}} \]

input
Integrate[(a + b*x^2)/(1 + x^2 + x^4)^2,x]
 
output
(x*(a + b - a*x^2 + 2*b*x^2))/(6*(1 + x^2 + x^4)) - (((-11*I + Sqrt[3])*a 
- 2*(-2*I + Sqrt[3])*b)*ArcTan[((-I + Sqrt[3])*x)/2])/(6*Sqrt[6 + (6*I)*Sq 
rt[3]]) - (((11*I + Sqrt[3])*a - 2*(2*I + Sqrt[3])*b)*ArcTan[((I + Sqrt[3] 
)*x)/2])/(6*Sqrt[6 - (6*I)*Sqrt[3]])
 
3.1.99.3 Rubi [A] (verified)

Time = 0.32 (sec) , antiderivative size = 128, normalized size of antiderivative = 1.08, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.389, Rules used = {1492, 1483, 1142, 25, 1083, 217, 1103}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {a+b x^2}{\left (x^4+x^2+1\right )^2} \, dx\)

\(\Big \downarrow \) 1492

\(\displaystyle \frac {1}{6} \int \frac {-\left ((a-2 b) x^2\right )+5 a-b}{x^4+x^2+1}dx+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

\(\Big \downarrow \) 1483

\(\displaystyle \frac {1}{6} \left (\frac {1}{2} \int \frac {5 a-b-3 (2 a-b) x}{x^2-x+1}dx+\frac {1}{2} \int \frac {5 a-b+3 (2 a-b) x}{x^2+x+1}dx\right )+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

\(\Big \downarrow \) 1142

\(\displaystyle \frac {1}{6} \left (\frac {1}{2} \left (\frac {1}{2} (4 a+b) \int \frac {1}{x^2-x+1}dx-\frac {3}{2} (2 a-b) \int -\frac {1-2 x}{x^2-x+1}dx\right )+\frac {1}{2} \left (\frac {1}{2} (4 a+b) \int \frac {1}{x^2+x+1}dx+\frac {3}{2} (2 a-b) \int \frac {2 x+1}{x^2+x+1}dx\right )\right )+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {1}{6} \left (\frac {1}{2} \left (\frac {1}{2} (4 a+b) \int \frac {1}{x^2-x+1}dx+\frac {3}{2} (2 a-b) \int \frac {1-2 x}{x^2-x+1}dx\right )+\frac {1}{2} \left (\frac {1}{2} (4 a+b) \int \frac {1}{x^2+x+1}dx+\frac {3}{2} (2 a-b) \int \frac {2 x+1}{x^2+x+1}dx\right )\right )+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

\(\Big \downarrow \) 1083

\(\displaystyle \frac {1}{6} \left (\frac {1}{2} \left (\frac {3}{2} (2 a-b) \int \frac {1-2 x}{x^2-x+1}dx-(4 a+b) \int \frac {1}{-(2 x-1)^2-3}d(2 x-1)\right )+\frac {1}{2} \left (\frac {3}{2} (2 a-b) \int \frac {2 x+1}{x^2+x+1}dx-(4 a+b) \int \frac {1}{-(2 x+1)^2-3}d(2 x+1)\right )\right )+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

\(\Big \downarrow \) 217

\(\displaystyle \frac {1}{6} \left (\frac {1}{2} \left (\frac {3}{2} (2 a-b) \int \frac {1-2 x}{x^2-x+1}dx+\frac {(4 a+b) \arctan \left (\frac {2 x-1}{\sqrt {3}}\right )}{\sqrt {3}}\right )+\frac {1}{2} \left (\frac {3}{2} (2 a-b) \int \frac {2 x+1}{x^2+x+1}dx+\frac {(4 a+b) \arctan \left (\frac {2 x+1}{\sqrt {3}}\right )}{\sqrt {3}}\right )\right )+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

\(\Big \downarrow \) 1103

\(\displaystyle \frac {1}{6} \left (\frac {1}{2} \left (\frac {(4 a+b) \arctan \left (\frac {2 x-1}{\sqrt {3}}\right )}{\sqrt {3}}-\frac {3}{2} (2 a-b) \log \left (x^2-x+1\right )\right )+\frac {1}{2} \left (\frac {(4 a+b) \arctan \left (\frac {2 x+1}{\sqrt {3}}\right )}{\sqrt {3}}+\frac {3}{2} (2 a-b) \log \left (x^2+x+1\right )\right )\right )+\frac {x \left (-\left (x^2 (a-2 b)\right )+a+b\right )}{6 \left (x^4+x^2+1\right )}\)

input
Int[(a + b*x^2)/(1 + x^2 + x^4)^2,x]
 
output
(x*(a + b - (a - 2*b)*x^2))/(6*(1 + x^2 + x^4)) + ((((4*a + b)*ArcTan[(-1 
+ 2*x)/Sqrt[3]])/Sqrt[3] - (3*(2*a - b)*Log[1 - x + x^2])/2)/2 + (((4*a + 
b)*ArcTan[(1 + 2*x)/Sqrt[3]])/Sqrt[3] + (3*(2*a - b)*Log[1 + x + x^2])/2)/ 
2)/6
 

3.1.99.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 217
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^( 
-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] & 
& (LtQ[a, 0] || LtQ[b, 0])
 

rule 1083
Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> Simp[-2   Subst[I 
nt[1/Simp[b^2 - 4*a*c - x^2, x], x], x, b + 2*c*x], x] /; FreeQ[{a, b, c}, 
x]
 

rule 1103
Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> S 
imp[d*(Log[RemoveContent[a + b*x + c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, 
e}, x] && EqQ[2*c*d - b*e, 0]
 

rule 1142
Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> S 
imp[(2*c*d - b*e)/(2*c)   Int[1/(a + b*x + c*x^2), x], x] + Simp[e/(2*c) 
Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x]
 

rule 1483
Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] : 
> With[{q = Rt[a/c, 2]}, With[{r = Rt[2*q - b/c, 2]}, Simp[1/(2*c*q*r)   In 
t[(d*r - (d - e*q)*x)/(q - r*x + x^2), x], x] + Simp[1/(2*c*q*r)   Int[(d*r 
 + (d - e*q)*x)/(q + r*x + x^2), x], x]]] /; FreeQ[{a, b, c, d, e}, x] && N 
eQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && NegQ[b^2 - 4*a*c]
 

rule 1492
Int[((d_) + (e_.)*(x_)^2)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symb 
ol] :> Simp[x*(a*b*e - d*(b^2 - 2*a*c) - c*(b*d - 2*a*e)*x^2)*((a + b*x^2 + 
 c*x^4)^(p + 1)/(2*a*(p + 1)*(b^2 - 4*a*c))), x] + Simp[1/(2*a*(p + 1)*(b^2 
 - 4*a*c))   Int[Simp[(2*p + 3)*d*b^2 - a*b*e - 2*a*c*d*(4*p + 5) + (4*p + 
7)*(d*b - 2*a*e)*c*x^2, x]*(a + b*x^2 + c*x^4)^(p + 1), x], x] /; FreeQ[{a, 
 b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && 
 LtQ[p, -1] && IntegerQ[2*p]
 
3.1.99.4 Maple [A] (verified)

Time = 0.16 (sec) , antiderivative size = 136, normalized size of antiderivative = 1.14

method result size
default \(-\frac {\left (\frac {a}{3}-\frac {2 b}{3}\right ) x -\frac {2 a}{3}+\frac {b}{3}}{4 \left (x^{2}-x +1\right )}-\frac {\left (6 a -3 b \right ) \ln \left (x^{2}-x +1\right )}{24}-\frac {\left (-2 a -\frac {b}{2}\right ) \sqrt {3}\, \arctan \left (\frac {\left (2 x -1\right ) \sqrt {3}}{3}\right )}{18}+\frac {\left (-\frac {a}{3}+\frac {2 b}{3}\right ) x -\frac {2 a}{3}+\frac {b}{3}}{4 x^{2}+4 x +4}+\frac {\left (6 a -3 b \right ) \ln \left (x^{2}+x +1\right )}{24}+\frac {\left (2 a +\frac {b}{2}\right ) \arctan \left (\frac {\left (1+2 x \right ) \sqrt {3}}{3}\right ) \sqrt {3}}{18}\) \(136\)
risch \(\frac {a \ln \left (124 a^{2} x^{2}-100 a b \,x^{2}+28 b^{2} x^{2}+124 a^{2} x -100 a b x +28 b^{2} x +124 a^{2}-100 a b +28 b^{2}\right )}{4}-\frac {b \ln \left (124 a^{2} x^{2}-100 a b \,x^{2}+28 b^{2} x^{2}+124 a^{2} x -100 a b x +28 b^{2} x +124 a^{2}-100 a b +28 b^{2}\right )}{8}+\frac {\left (-\frac {a}{6}+\frac {b}{3}\right ) x^{3}+\left (\frac {a}{6}+\frac {b}{6}\right ) x}{x^{4}+x^{2}+1}+\frac {\sqrt {3}\, a \arctan \left (\frac {62 \sqrt {3}\, a^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {50 \sqrt {3}\, a x b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {14 \sqrt {3}\, b^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {31 \sqrt {3}\, a^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {25 \sqrt {3}\, a b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {7 \sqrt {3}\, b^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}\right )}{9}+\frac {\sqrt {3}\, b \arctan \left (\frac {62 \sqrt {3}\, a^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {50 \sqrt {3}\, a x b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {14 \sqrt {3}\, b^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {31 \sqrt {3}\, a^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {25 \sqrt {3}\, a b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {7 \sqrt {3}\, b^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}\right )}{36}-\frac {a \ln \left (124 a^{2} x^{2}-100 a b \,x^{2}+28 b^{2} x^{2}-124 a^{2} x +100 a b x -28 b^{2} x +124 a^{2}-100 a b +28 b^{2}\right )}{4}+\frac {b \ln \left (124 a^{2} x^{2}-100 a b \,x^{2}+28 b^{2} x^{2}-124 a^{2} x +100 a b x -28 b^{2} x +124 a^{2}-100 a b +28 b^{2}\right )}{8}+\frac {\sqrt {3}\, a \arctan \left (\frac {62 \sqrt {3}\, a^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {50 \sqrt {3}\, a x b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {14 \sqrt {3}\, b^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {31 \sqrt {3}\, a^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {25 \sqrt {3}\, a b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {7 \sqrt {3}\, b^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}\right )}{9}+\frac {\sqrt {3}\, b \arctan \left (\frac {62 \sqrt {3}\, a^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {50 \sqrt {3}\, a x b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {14 \sqrt {3}\, b^{2} x}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {31 \sqrt {3}\, a^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}-\frac {25 \sqrt {3}\, a b}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}+\frac {7 \sqrt {3}\, b^{2}}{3 \left (31 a^{2}-25 a b +7 b^{2}\right )}\right )}{36}\) \(906\)

input
int((b*x^2+a)/(x^4+x^2+1)^2,x,method=_RETURNVERBOSE)
 
output
-1/4*((1/3*a-2/3*b)*x-2/3*a+1/3*b)/(x^2-x+1)-1/24*(6*a-3*b)*ln(x^2-x+1)-1/ 
18*(-2*a-1/2*b)*3^(1/2)*arctan(1/3*(2*x-1)*3^(1/2))+1/4*((-1/3*a+2/3*b)*x- 
2/3*a+1/3*b)/(x^2+x+1)+1/24*(6*a-3*b)*ln(x^2+x+1)+1/18*(2*a+1/2*b)*arctan( 
1/3*(1+2*x)*3^(1/2))*3^(1/2)
 
3.1.99.5 Fricas [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 185, normalized size of antiderivative = 1.55 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=-\frac {12 \, {\left (a - 2 \, b\right )} x^{3} - 2 \, \sqrt {3} {\left ({\left (4 \, a + b\right )} x^{4} + {\left (4 \, a + b\right )} x^{2} + 4 \, a + b\right )} \arctan \left (\frac {1}{3} \, \sqrt {3} {\left (2 \, x + 1\right )}\right ) - 2 \, \sqrt {3} {\left ({\left (4 \, a + b\right )} x^{4} + {\left (4 \, a + b\right )} x^{2} + 4 \, a + b\right )} \arctan \left (\frac {1}{3} \, \sqrt {3} {\left (2 \, x - 1\right )}\right ) - 12 \, {\left (a + b\right )} x - 9 \, {\left ({\left (2 \, a - b\right )} x^{4} + {\left (2 \, a - b\right )} x^{2} + 2 \, a - b\right )} \log \left (x^{2} + x + 1\right ) + 9 \, {\left ({\left (2 \, a - b\right )} x^{4} + {\left (2 \, a - b\right )} x^{2} + 2 \, a - b\right )} \log \left (x^{2} - x + 1\right )}{72 \, {\left (x^{4} + x^{2} + 1\right )}} \]

input
integrate((b*x^2+a)/(x^4+x^2+1)^2,x, algorithm="fricas")
 
output
-1/72*(12*(a - 2*b)*x^3 - 2*sqrt(3)*((4*a + b)*x^4 + (4*a + b)*x^2 + 4*a + 
 b)*arctan(1/3*sqrt(3)*(2*x + 1)) - 2*sqrt(3)*((4*a + b)*x^4 + (4*a + b)*x 
^2 + 4*a + b)*arctan(1/3*sqrt(3)*(2*x - 1)) - 12*(a + b)*x - 9*((2*a - b)* 
x^4 + (2*a - b)*x^2 + 2*a - b)*log(x^2 + x + 1) + 9*((2*a - b)*x^4 + (2*a 
- b)*x^2 + 2*a - b)*log(x^2 - x + 1))/(x^4 + x^2 + 1)
 
3.1.99.6 Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.93 (sec) , antiderivative size = 874, normalized size of antiderivative = 7.34 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=\text {Too large to display} \]

input
integrate((b*x**2+a)/(x**4+x**2+1)**2,x)
 
output
(x**3*(-a + 2*b) + x*(a + b))/(6*x**4 + 6*x**2 + 6) + (-a/4 + b/8 - sqrt(3 
)*I*(4*a + b)/72)*log(x + (76*a**3*(-a/4 + b/8 - sqrt(3)*I*(4*a + b)/72) + 
 948*a**2*b*(-a/4 + b/8 - sqrt(3)*I*(4*a + b)/72) - 816*a*b**2*(-a/4 + b/8 
 - sqrt(3)*I*(4*a + b)/72) + 12096*a*(-a/4 + b/8 - sqrt(3)*I*(4*a + b)/72) 
**3 + 148*b**3*(-a/4 + b/8 - sqrt(3)*I*(4*a + b)/72) - 8640*b*(-a/4 + b/8 
- sqrt(3)*I*(4*a + b)/72)**3)/(248*a**4 - 262*a**3*b + 75*a**2*b**2 + 11*a 
*b**3 - 7*b**4)) + (-a/4 + b/8 + sqrt(3)*I*(4*a + b)/72)*log(x + (76*a**3* 
(-a/4 + b/8 + sqrt(3)*I*(4*a + b)/72) + 948*a**2*b*(-a/4 + b/8 + sqrt(3)*I 
*(4*a + b)/72) - 816*a*b**2*(-a/4 + b/8 + sqrt(3)*I*(4*a + b)/72) + 12096* 
a*(-a/4 + b/8 + sqrt(3)*I*(4*a + b)/72)**3 + 148*b**3*(-a/4 + b/8 + sqrt(3 
)*I*(4*a + b)/72) - 8640*b*(-a/4 + b/8 + sqrt(3)*I*(4*a + b)/72)**3)/(248* 
a**4 - 262*a**3*b + 75*a**2*b**2 + 11*a*b**3 - 7*b**4)) + (a/4 - b/8 - sqr 
t(3)*I*(4*a + b)/72)*log(x + (76*a**3*(a/4 - b/8 - sqrt(3)*I*(4*a + b)/72) 
 + 948*a**2*b*(a/4 - b/8 - sqrt(3)*I*(4*a + b)/72) - 816*a*b**2*(a/4 - b/8 
 - sqrt(3)*I*(4*a + b)/72) + 12096*a*(a/4 - b/8 - sqrt(3)*I*(4*a + b)/72)* 
*3 + 148*b**3*(a/4 - b/8 - sqrt(3)*I*(4*a + b)/72) - 8640*b*(a/4 - b/8 - s 
qrt(3)*I*(4*a + b)/72)**3)/(248*a**4 - 262*a**3*b + 75*a**2*b**2 + 11*a*b* 
*3 - 7*b**4)) + (a/4 - b/8 + sqrt(3)*I*(4*a + b)/72)*log(x + (76*a**3*(a/4 
 - b/8 + sqrt(3)*I*(4*a + b)/72) + 948*a**2*b*(a/4 - b/8 + sqrt(3)*I*(4*a 
+ b)/72) - 816*a*b**2*(a/4 - b/8 + sqrt(3)*I*(4*a + b)/72) + 12096*a*(a...
 
3.1.99.7 Maxima [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 105, normalized size of antiderivative = 0.88 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=\frac {1}{36} \, \sqrt {3} {\left (4 \, a + b\right )} \arctan \left (\frac {1}{3} \, \sqrt {3} {\left (2 \, x + 1\right )}\right ) + \frac {1}{36} \, \sqrt {3} {\left (4 \, a + b\right )} \arctan \left (\frac {1}{3} \, \sqrt {3} {\left (2 \, x - 1\right )}\right ) + \frac {1}{8} \, {\left (2 \, a - b\right )} \log \left (x^{2} + x + 1\right ) - \frac {1}{8} \, {\left (2 \, a - b\right )} \log \left (x^{2} - x + 1\right ) - \frac {{\left (a - 2 \, b\right )} x^{3} - {\left (a + b\right )} x}{6 \, {\left (x^{4} + x^{2} + 1\right )}} \]

input
integrate((b*x^2+a)/(x^4+x^2+1)^2,x, algorithm="maxima")
 
output
1/36*sqrt(3)*(4*a + b)*arctan(1/3*sqrt(3)*(2*x + 1)) + 1/36*sqrt(3)*(4*a + 
 b)*arctan(1/3*sqrt(3)*(2*x - 1)) + 1/8*(2*a - b)*log(x^2 + x + 1) - 1/8*( 
2*a - b)*log(x^2 - x + 1) - 1/6*((a - 2*b)*x^3 - (a + b)*x)/(x^4 + x^2 + 1 
)
 
3.1.99.8 Giac [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 109, normalized size of antiderivative = 0.92 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=\frac {1}{36} \, \sqrt {3} {\left (4 \, a + b\right )} \arctan \left (\frac {1}{3} \, \sqrt {3} {\left (2 \, x + 1\right )}\right ) + \frac {1}{36} \, \sqrt {3} {\left (4 \, a + b\right )} \arctan \left (\frac {1}{3} \, \sqrt {3} {\left (2 \, x - 1\right )}\right ) + \frac {1}{8} \, {\left (2 \, a - b\right )} \log \left (x^{2} + x + 1\right ) - \frac {1}{8} \, {\left (2 \, a - b\right )} \log \left (x^{2} - x + 1\right ) - \frac {a x^{3} - 2 \, b x^{3} - a x - b x}{6 \, {\left (x^{4} + x^{2} + 1\right )}} \]

input
integrate((b*x^2+a)/(x^4+x^2+1)^2,x, algorithm="giac")
 
output
1/36*sqrt(3)*(4*a + b)*arctan(1/3*sqrt(3)*(2*x + 1)) + 1/36*sqrt(3)*(4*a + 
 b)*arctan(1/3*sqrt(3)*(2*x - 1)) + 1/8*(2*a - b)*log(x^2 + x + 1) - 1/8*( 
2*a - b)*log(x^2 - x + 1) - 1/6*(a*x^3 - 2*b*x^3 - a*x - b*x)/(x^4 + x^2 + 
 1)
 
3.1.99.9 Mupad [B] (verification not implemented)

Time = 13.59 (sec) , antiderivative size = 897, normalized size of antiderivative = 7.54 \[ \int \frac {a+b x^2}{\left (1+x^2+x^4\right )^2} \, dx=\text {Too large to display} \]

input
int((a + b*x^2)/(x^2 + x^4 + 1)^2,x)
 
output
atan((((2*b - 10*a + 24*x*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/ 
72))*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) - x*((59*a^2)/18 
- (19*a*b)/9 + b^2/9))*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) 
*1i + ((10*a - 2*b + 24*x*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/ 
72))*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) - x*((59*a^2)/18 
- (19*a*b)/9 + b^2/9))*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) 
*1i)/((19*a*b^2)/36 - (29*a^2*b)/36 + (31*a^3)/108 - (7*b^3)/54 + ((2*b - 
10*a + 24*x*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72))*(b/8 - a/ 
4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) - x*((59*a^2)/18 - (19*a*b)/9 + 
 b^2/9))*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) - ((10*a - 2* 
b + 24*x*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72))*(b/8 - a/4 + 
 (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72) - x*((59*a^2)/18 - (19*a*b)/9 + b^ 
2/9))*(b/8 - a/4 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72)))*((a*1i)/2 - (b 
*1i)/4 + (3^(1/2)*a)/9 + (3^(1/2)*b)/36) + atan((((2*b - 10*a + 24*x*(a/4 
- b/8 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72))*(a/4 - b/8 + (3^(1/2)*a*1i 
)/18 + (3^(1/2)*b*1i)/72) - x*((59*a^2)/18 - (19*a*b)/9 + b^2/9))*(a/4 - b 
/8 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72)*1i + ((10*a - 2*b + 24*x*(a/4 
- b/8 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72))*(a/4 - b/8 + (3^(1/2)*a*1i 
)/18 + (3^(1/2)*b*1i)/72) - x*((59*a^2)/18 - (19*a*b)/9 + b^2/9))*(a/4 - b 
/8 + (3^(1/2)*a*1i)/18 + (3^(1/2)*b*1i)/72)*1i)/((19*a*b^2)/36 - (29*a^...